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Abstract 

The capillary-induced structural instability of an elastic circular tube partially filled by a liquid is 
studied by combining theoretical analysis and molecular dynamics simulations. The analysis shows 
that, associated with the instability, there is a well-defined length scale (elasto-capillary length), 
which exhibits a scaling relationship with the characteristic length of the tube, regardless of the 
interaction details. We validate this scaling relationship for a carbon nanotube partially filled by 
liquid iron. The capillary-induced structural transformation could have potential applications for 
nano-devices. 
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A wide range of phenomena in nature, which span from everyday observations to many 
bio-related processes, are a result of capillary forces. Common examples include the shape 
of liquid droplets, the imbibition of a sponge the clumping of wet hair into bundles 2] 

and coalescence of paintbrush fibers j3|, the standing of aquatic insects on water [4l|, and 

I, 

lung airway closure etc |5|]. Recently it was found that capillary forces can induce structural 
deformations or instability in an elastic system, in cases where these surface forces are 
comparable in magnitude to bulk elastic forces. More importantly, there is a well-defined 
length scale (elasto-capillary length Lec) which underlies such structural instabilities 6|, 



which usually reveals an intrinsic scaling relation with the characteristic dimensions 



a. 
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The elasto-capillarity deformation or instability has been also found to have many 



interesting applications 
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Liquid infiltrated elastic tubes are common in bio-systems and everyday life. When the 
capillary forces are comparable to the bending stiffness of the tube, the tube can be deformed 
or become unstable at some critical filling fraction. Whether the instability can be associated 
with a new type of scaling relationship has thus far remained relatively unexplored. In this 
work, molecular dynamics (MD) simulations of this phenomenon are coupled with theoretical 
analysis to investigate this issue for an elastic cylindrical tube partially filled by a liquid. 
Our results show that the instability or deformation of a flexible tube leads to a well- 
defined scaling law regardless of the interaction details. Our MD simulations on liquid iron 
encapsulated by carbon nanotubes quantitatively demonstrate the scaling law. 

We consider an elastic circular tube with length of Lt and radius of Rq partially filled 
with an incompressible Newtonian liquid of length L;, where Rq refers the radius of the 
"dry" tube without any compression. The upper panel of Fig. 1 illustrates the model used 
in the current work. The energy of the system can be written formally as, 

E = E, + E, + Ei (1) 

where Ef,, E^ and Ei are the elastic bending strain energy, elastic compression strain energy, 
and interfacial free energy (capillary energy), respectively. In the continuum elasticity pic- 



ture 
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131 ] , the bending strain energy is = ^ / ^dl, and the compression strain energy 



is = M: fii^y'dl, where / denotes the curvliuea.^ integral along the pernneter of the 

cross section and Pq = 2ttRq is the perimeter of the original tube cross section. The bending 
stiffness B = Yh^ / 12(1 — u"^) and the compression stiffness C = Yh/{l — v'^) are constants re- 
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lated to Young's modulus Y, the Poisson ratio z/, and the tube thickness h. The local radius 
of curvature is denoted as p. The interfacial free energy is, Ei = Ai^'yiy+Asi'jsi+Asv'ysv+'^PT, 
where Aiy, Asi and Ag^ are the areas of liquid- vapor, solid-liquid and solid- vapor interfaces, 
respectively, P is the perimeter of the cross section and Ag^ = 2PLt — Agi, here we have 
included the solid-vapor surface outside the cylinder that is not covered by liquid. The sur- 
face and interface free energies are denoted 7;^, and 7s^,, for the hquid- vapor, solid-liquid 
and solid-vapor interfaces, respectively. The specific free energy of the three-phase contact 
line (i.e., the line tension) is denoted as r.[l^ The 2Pr term usually can be neglected if the 
dimension of the system is not too small. 

Since only near-critical sizes (where the circular shape is near its stability limit) are 
considered, the non-circular shape is approximated by that of an ellipse. For an ellipse 
with long axis Ra and short axis Rb-, the shape can be described well by introducing two 
independent variables, R and tu, which are defined as i? = \/RaRb and u = Ra/ Rb- When 
u!=l, these variables describe a perfect circle. By introducing two integrals, fi,2{^) = 
^ f{u cos"^ t + u^^ sin^ t)°'^''^dt, with ai = —2.5 and 0:2 = 0.5, we can write / jsdl = ^fi{u) 
and § dl = P = 27ri?/2(a;).[15| Finally we have 

Eb + E, = ^Mco) + CtcRoL.H^Uco))' - 2{^Uu)f + ^/^(c.)]. (2) 

li IIQ IhQ IhQ 

In order to calculate the capillary energy, the area of the interfaces should be written 
in explicit forms. However, the general calculation of the liquid-vapor interface, which is a 
meniscus, for arbitrary tube shape is beyond the scope of the present treatment and will 
be discussed in a subsequent paper [l6|. To elucidate the important physics as well as 
to simplify the mathematical treatment, we consider the cases where the contact angles are 
close to 90°. Within this limitation, all the interface areas can be defined without ambiguity. 
It should be noted that in the simulations discussed below, the solid-liquid contact angles 
are in fact approximately 90°. 

Thermal fluctuations of the tubes are neglected as they exist on a length scale compa- 
rable to the persistence length(/p = YI/ksT, and Y: Young's modulus, I: the area moment 
of intertia. ), which usually is much longer than the characteristic length of system. For 
example, in single-walled carbon nanotubes, the persistence length is on the order of 45 fim 
6j. With the above assumptions as well as the assumption of that the liquid is incompress- 
ible, we can write Ai^ = 27iR , Asi = PLi{-^) and A^y = 2PLt — A^i. So, the interfacial free 
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energy is, 

E, = 27tR\, + PLi^ilsi - 7s.) + 2PLasv + 2Pt (3) 

In Eq. 2, the elastic energy always tends to maintain the tube in a circular shape. If 
the capillary energy approaches or exceeds the magnitude of the elastic energy, the circular 
shape will be unstable and the system might seek a lower energy state. Physically, the 
deformation of the tube reduces the interfacial energy at the cost of increasing the elastic 
energy. For 9 < 90°, (7^; — 7^^,) < 0, the system could deform to increase the area of the 
solid-liquid interface. Even for 6 > 90°, (7^/ — 7si.) > 0, the tube still has a chance to deform 
if the first term exceeds the second term in Eq. 3, in which case the area of the liquid- vapor 
interface decreases. The MD simulations in the current work focus on the latter case. The 
critical transition point can be located by the condition ^^\ui=i = 0. And the radius of 
tube at the critical point (Re) is determined by U=i = 0. Remembering /i,2(^^)U=i = 1, 
/i,2('-^)U=i — /i'2('-^)U=i — I' above conditions produce the critical relation. 



„ 2|7,„-(7.,-7„,).ii|i].fl^ 

^' = w ' 

where is the critical length of the tube. For elastic tubes, if the compression stiffness (C) 
is much larger than the bending stiffness B, which is the the case we are interested in, one 
can neglect the different between Rq and Rc, 

r._ ^hv-hsl-%.)-^]-Rl 

~ w 

We can define an effective interface free energy 7* = [7;^ — (7^/ — 7s„) ■ ^], where the 



factor in 7* reflects an inherent geometric relation among the three interfaces. Note 



Ro 



that if 27* is divided L*, according to Eq. 5, a generalized pressure p* = ^ = ^ can be 



obtained, which is the same form as reported in Ref. jl5| for CNTs under pressure. In fact, 
this capillary force can be regarded as a kind of negative internal pressure. From Eq. 5, a 
generalized elasto-capillary scaling relation can be derived. 
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Similar to previous studies, we can define an elasto-capillary length {Lec = ["^] 2) 

for this system, which gives the typical effective curvature induced by capillarity on the 
tube. Different from previous studies, the elasto-capillary length is not simply defined in 



terms of interfacial free energies, but rather an effective (average) interface free energy (7*). 
The difference stems from the fact that in the filled liquid tube system studied here, the 
area of the solid-liquid (solid-vapor) and liquid-vapor interfaces changes in different ways 
(amounts) as the tube is deformed. The former scales as R, while the later scales as R^. It 
implies that the average capillary forces depend on the system size itself (-^^ in 7* reflects the 
dependence). However in previous studies, either the three interfaces change simultaneously 
by the same amount (liquid drop on thin films 

n n n 

(e.g., slender rods immersed in liquid) [2, 13|, 16|, 
independent on the size of system. 

We can also define a characteristic length of the system at the critical point, Lc = 

Area of cross section 2ttrI z^— /2_r3x1 1 • 1 • , • j • j.- r 
— = voTT^oT*-) 2 , which IS a purely geometric description of 



or only one interface changes area 
9|, and thus the average capillary force is 



^Surface area of tube ^J2-kRoLi ^""^3L_ 
the tube. Eq. 6 becomes a generic scaling relationship 



EC 
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Lq defines a typical length scale at which the elastic energy and surface energy are compa- 
rable. Different from previous studies for slender rods immersed in liquid, j^,^, [g], or for 
liquid drops on thin films, 7, ^ the characteristic length in the present case is not simply 
a spatial dimension of system. The reason is that the elastic energy depends on both Rq 
and Lf, where Rq and Lt are functions of the radial bulk modulus and total elastic energy 
respectively. Therefore, the elasto-capillary length and the scaling relationship presented 
here represent a new type of instability process. 

To test the theoretical analysis, we simulate liquid iron encapsulated by a single-walled 
carbon nanotube (SWCNT) using MD simulation. All the simulations made use of the 
LAMMPS (large-scale atomic/molecular massively parallel simulator) code |17i]. The inter- 
action between carbon atoms is described by the second-generation reactive empirical bond 
order (REB02) potential. [3] The niany-body potential developed by Mendelev et al is used 
to describe the Fe-Fe interactions. [19] The weak interaction between C and Fe is modeled 



by a truncated Lennard-Jones(LJ) potential {20], |2l|. The interaction parameters a and e 
are obtained by fitting to ab-initio results, {22] a = 2.05A ande = 0.09463eV. The ab-initio 
results are also used by other authors to fit the C-Fe interaction based on Johnson pair 



form potential. 23|] This fitted potential give a reasonable contact angle (~ 107" at 2500K) 



compared to the experimental value (~ 125° at 923K) for liquid Fe in SWCNT 2^ . 
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In the present study, periodic boundary conditions in the axial direction and free bound- 
ary conditions in the radial directions are adopted. The systems are simulated at T=2500 
K and zero pressure. The temperature is much higher than the melting point of Fe , 19| so 
that Fe is in a well defined liquid state. A typical run lasts 3ns for each case. The tubes 
used in this work are of the zigzag type, including (32,0), (40,0), (46,0) and (50,0), with the 
typical length in the range of 30A to 250A. To locate the critical length for each sample, 
the length of the tube is decreased until the circular shape of the CNT becomes unstable. 

To quantitatively compare the model prediction and MD results, the value of all the 
variables on the right hand side of Eq. 5 are needed. -y/„ for liquid iron at 2500K is 



calculated separately using the method reported in Refs. 25|,|26||, the value obtained is about 



0.82±0. 02 J/m^. Both {jsi—lsv) and r can be calculated from the extended Young's equation 



14| if the contact angle as a function of the radius of the CNT is known. Contact angles 6 



are measured during MD simulations for a circular shape. The determination of the contact 
angle follows the procedure described in Refs. 27, '2^- We obtain 7^/ — 7st,=0.24±0. 02 J/m^ 
and r = 0.4 ± 0.3 x 10^^° J/m. The elastic constant B of the CNT for the REB02 potential 



is adopted from Ref. 



29|, B=lAeV. 



The measurement of the CNT shape is based on the calculation of the cross sectional area 
S. An order parameter thus can be defined as A =< S >t /So, where So is the initial cross 
sectional area of the CNT and < ... >t denotes the time average. With this definition A is 
one for the circular shape while it is smaller than one for a non-circular shape, and decreases 
continuously with increasing deviation from a circular shape. Due to thermal fluctuations, 
the shape of the CNT varies with time. We define the critical length for which A is smaller 
than one within the statistical uncertainty. 

The lower panel of Fig. 1 shows the time-averaged cross section for a (40,0) CNT filled 
with I6.85A liquid for a few selected Lt. With decreasing Lt, the shape of the cross section 
changes from a circle to an ellipse. Fig. 2 shows the evolution of A vs the length of the CNT 
{Lt) for a (40,0) tube, where the different symbols denote systems with different lengths 
of liquid. We find that the order parameter can clearly distinguish the CNT shape change 
from circle to ellipse. From Fig. 2, one can see that, for longer Lf, A remains at an average 
value of one (corresponding to a circular shape) within the error bars, while A deviates from 
one for shorter Lt. The critical length can be located from this figure. It should be pointed 
out that the critical length is not dependent on the specific definition of order parameter. 
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;her order parameter based on the Fourier spectrum of a shape fluctuation function 



31| was also tested and the results do not show signiflcant differences. For all studied 



samples, when Lf is smaller than a certain critical length, the shape of the cross section will 
change from circle to ellipse. The shape changes are purely elastic deformations and the 
system can relax reversibly by removing the liquid. 

Fig. 3 shows the critical length of the tube vs the liquid length Li for all CNTs studied, 
where the solid line is plotted based on Eq. 5 using the parameter values given above. One 
can see that the MD results are in excellent agreement with the theoretical predictions. The 
MD results are consistent with a linear relationship between Lc and Leg as expected from 
the model analysis (see Fig. 4). A best flt produces a slope of 3.034±0.014, which is in 
excellent agreement with the theoretical prediction of ^/Stt. 

The critical length (Eq. 5) affords a few potential applications in nano-science or in 
designing nano-devices. One possible application is the measurement of the elastic bending 
stiffness (B) of the tube or the contact angle. For this purpose, we only need knowledge of 
as a function of Li by fllling liquid into the tubes. If the and contact angle are known, B 
can be obtained. On the other hand, if B is known, the contact angle can be estimated. The 
capillary induced structural transformation also provides an ideal mechanism of conversion 
of capillary energy to elastic energy. If the values of the interface energies in Eq. 5 are 
sensitive to the temperature, or to other environment variables, an environment-controlled 
nano-machine could be achieved. 

In summary, the capillary-induced instability of elastic circular tubes has been studied 
by combining theoretical analysis and MD simulation. The instability depends on two key 
length-scales, the elasto-capillary length Lec which represents the typical radius of curvature 
produced by capillary forces on an elastic tube and the characteristic geometric length Lc, 
which deflnes a typical reduced size of the tube. At the critical point, Lc is found to be 
proportional to Lec, above which a tube preserves a circular shape. The MD results are 
found to be in excellent agreement with the theoretical prediction. The current results can 
be useful to the design of nano-devices. 

We thank Prof. B. B. Laird for his helpful discussion. This research is supported by the 
National Science Foundation of China, the special funds for major state basic research and 
CAS projects. 
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FIG. 1. (Color online) Upper panel: Sketch of the model used in present work. The elastic tube is 
characterized by length [Lt) and radius {Rq). The encapsulated liquid is described by the length 
(L;), the height of cap (r/) and the contact angle {6). Lower panel: Time averaged cross section for 
various lengths of the CNT Lt. All results are for (40,0) CNTs with hquid length L/ = 16.851. (a) 
Lt = 96.921, (b) Lt = 84.281, (c) Lt = 59.00l, (d) Lt = 50.581, (e) Lt = 46.361, (f) Lt = 42.151. 
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FIG. 3. The critical length of the tube as a function of the length of liquid filled region. Symbols: 
the molecular dynamics results, solid lines: the analysis results of Eq. 5, dashed lines: The tube is 
fully filled with liquid, below this line the model is not valid. 
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FIG. 4. The characteristic geometry length as a function of elasto-capiUary length. Symbols: the 
molecular dynamics results, lines: analysis results of Eq. 7. 
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